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ABSTRACT 

Recently, it was suggested that a photospheric component that resuhs from the internal dissipation occurring 
in the optically thick inner parts of relativistic outflows may be present in the prompt 7/X-ray emission of 
gamma-ray bursts or X-ray flashes. We explore high-energy neutrino emission in this dissipative photosphere 
model, assuming that the composition of the outflow is baryon-dominated. We find that neutrino emission 
from proton-proton collision process forms an interesting signature in the neutrino spectra. Under favorable 
conditions for the shock dissipation site, these low-energy neutrinos could be detected by km^ detectors, such 
as Icecube. Higher energies (> 10 TeV) neutrino emission from proton-proton collision and photo-pion pro- 
duction processes could be significantly suppressed for dissipation at relatively small radii, due to efficient 
Bethe-Heitler cooling of protons and/or radiative cooling of the secondary mesons in the photosphere radia- 
tion. As the dissipation shocks continue further out, high energy neutrinos from photo-pion production process 
becomes dominant. 

Subject headings: gamma rays: bursts; elementary particles 



1. INTRODUCTION 

Although it has been generally accepted that the prompt 
gamma-ray emission of gamma-ray bursts (GRBs) results 
from internal dissipation, likely internal shocks, of a rel- 
ativistic outflow (e.g. Paczyhski & Xu 1994; Rees & 
Meszaros 1994), the dissipation site and the radiation mech- 
anism for the gamma-ray emission are still largely unknown. 
Synchrotron and/or inverse-Compton scattering emission by 
shock-accelerated electrons in the optically-thin region has 
been proposed as an efficient mechanism for the gamma-ray 
emission. However, this model does not satisfactorily account 
for a few observational facts, such as the low-energy spectral 
slops that are steeper than synchrotron lower energy spectral 
indices (Preece et al. 2000; Lloyd et al. 2000), the cluster- 
ing of peak energies, the correlation between the burst's peak 
energy and luminosity (Amati et al. 2000). It becomes rec- 
ognized that an additional thermal component may play a key 
role and could solve these problems (e.g. Pe'er et al. 2006; 
Ryde et al. 2006). It has also been pointed out that a hy- 
brid model with both a thermal and non-thermal component 
can describe the spectrum equally well as the Band function 
model (Band 1993), but the former has a more physical mean- 
ing (Ryde 2005). Recently, it was suggested that a strong 
quasi-thermal component could result from the internal dissi- 
pation occurring in the optically thick inner parts of relativis- 
tic outflows (Rees & Meszaros 2005; Pe'er, Meszaros & Rees 
2006; Thompson et al. 2007). Sub-photospheric shock dissi- 
pation can increase the radiative efficiency of the outflow, sig- 
nificantly boosting the original thermal photospheric compo- 
nent so that it may well dominate the nonthermal component 
from optically-thin shocks occurring outside the photosphere. 

Neutrino emission from gamma-ray bursts has been pre- 
dicted at different stages of the relativistic outflow, such as the 
precursor phase (e.g. Bahcall & Meszaros 2000; Meszaros & 
Waxman 2001; Razzaque et al. 2003a,b; Razzaque et al. 
2004; Ando & Beacom 2005; Horiuchi & Ando 2008; Ko- 
ers & Wijers 2008), the prompt emission phase (e.g. Waxman 
& Bahcall 1997; Dermer & Atoyan 2003; Guetta et al. 2004; 
Murase & Nagataki 2006; Gupta & Zhang 2007; Murase et 
al. 2006) and afterglow phase (e.g. Waxman & Bahcall 2001; 



Dai & Lu 2001; Dermer 2002; Li et al. 2002; Murase & Na- 
gataki 2006; Murase 2007; Dermer 2007). Based on the bro- 
ken power-law approximation for the spectrum of the prompt 
emission, presumably from optically-thin internal shocks, a 
burst of PeV neutrinos, produced by photomeson production, 
was predicted to accompany the prompt gamma-ray emission 
if protons are present and also accelerated in the shocks (Wax- 
man & Bahcall 1997). The neutrino emission from proton- 
proton ipp) collisions was generally thought to be negligi- 
ble due to lower collision opacity for optically-thin internal 
shocks. However, as we show below, if some part of the 
prompt emission arises from internal shocks occurring in the 
optically-thick inner part of the outflow, as indicated by the 
thermal emission, a lower energy (< lOTeV) neutrino com- 
ponent may appear as a result of pp collisions. 

2. THE DISSIPATIVE PHOTOSPHERE MODEL 

The photosphere models have been widely discussed in re- 
lation to the prompt emission of GRBs (e.g. Thompson 1994; 
Ghiseflini & Celotti 1999; Rees & Meszaros 2005; Thomp- 
son et al. 2007; loka et al. 2007). The potential advantage 
of photosphere models is that the peak energy can be stabi- 
lized, which is identified as the thermal or Comptonization 
thermal peak (see loka et al. 2007 and references therein). 
The photosphere radiation may also produce a large num- 
ber of electron-positron pairs, which may lead to a pair pho- 
tosphere beyond the baryon-related photosphere (e.g. Rees 
& Meszaros 2005), and may also enhance the radiative effi- 
ciency (loka et al. 2007). On the other hand, it is also sug- 
gested that the number of pairs produced does not exceed the 
baryon related electrons by a factor larger than a few (Pe'er 
et al. 2006). For simplicity, we here only consider the dissi- 
pation below the baryon-related photosphere, which is more 
favorable for pp neutrino production. 

Foflowing Rees & Meszaros (2005) and Pe'er et al. (2006), 
we assume that during the early stage of the prompt emission, 
internal shocks of the outflow occur at radii below the bary- 
onic photosphere. Initially, the internal energy is released at 
the base of the outflow, ro ^ arg = laGM /c^, where a > 1 
and rg is the Schwarzschild radius of a central object of mass 
M. The internal energy is then converted to the kinetic energy 
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of the flow, whose bulk Lorentz factor grows as 7 ~ r up to a 
saturation radius at ~ ror], where = Lq /(Mc^) is the initial 
dimensionless entropy, Lq and M are the total energy and mass 
outflow rates. Above the saturation radius, the observer-frame 
photospherical luminosity decreases as Ly(r) = Loir / r^)'^^^ 
and the greater part of energy is in kinetic form, ~ Lq. If 
the dissipation is maintained all the way to the photosphere, it 
will lead to an effective luminosity Lj ~ e^Lo and a tempera- 
ture (Rees & Meszaros 2005) 



r^rfLl/^keV, (1) 



where 64 is the dissipation efficiency, Tq is the initial temper- 
ature of the fireball outflow and F is the bulk Lorentz fac- 
tor of the outflow. The internal shock occurs at /? ~ ITr^ = 
6 X lOi°ar^(M/lOM0)cm, for which the optical depfli to 
Thomson scattering by the baryon-related electrons is tj = 
aTLk/(47rRT^mpC^) = l20Lk,52Ru^2^ , where aj is the Thom- 
son scattering cross section and is the kinetic energy lu- 
minosity. The photosphere is further out, at radius Rpt, = 
1.2 X lO^^Lir saFj^cm. A detailed calculation taking into ac- 
count of the electron/positron cooling and the Comptonization 
effect leads to a quasi-thermal emission which peaks at energy 
^ 300-500keV for dissipation at Thomson optical depth of 
tt ~ 10-100 (Pe'er et al. 2006). This temperature is con- 
sistent with the observed peak energies of prompt gamma-ray 
emission of a majority of GRBs. 

Assuming that a fraction of 0.1 of the shock inter- 
nal energy is converted into magnetic fields, we have a mag- 
netic field B' = 2.5 X lO^e^^jL^^jj/^Tjrj'G. Protons accel- 
erated by internal shocks are assumed to have a spectrum 
dn/dsp ~ EpP with ~ 2, as often assumed for non-relativistic 
or nuldly-relativistic shock acceleration. The maximum pro- 
ton energy is set by comparing the acceleration time scale 

-12^.^^V-|/fL-'/2'/?„r2Swith 



C = ae'JieB'c) = 4.4 x lO-i^a(j^)e 



the energy-loss time scales. The synchrotron loss time is 



Assuming that the sub-photosphere emission at the dissipa- 
tion site peaks at = 300keV with a thermal-like spectrum, 
the number density of photons in the comoving frame is niy = 
L^/i4wR^T^cs'^) = 5 X lO^^L^^siRf^r^^e^/SOOkeYT^cm-^. 
The pj cooling time is approximately f^^ = 1 /{(jp^n'^cKp~f) = 
10~"^L~'5j/?j J 3ookev ) whcrc Kp^ ~ 0.2 is the inelasticity 
and cTp^ = 5 X lO'^'^cm^ is the peak cross section at the A 
resonance. By comparison with the synchrotron loss time, it 
is found that the most effective coohng mechanism for pro- 
tons is the P'y process for protons with energies above the p'y 
threshold, but below ~ lO^GeV. Equating t'^^^ = tp^, we obtain 
the maximum proton energy in the shock comoving frame 

)GeV. (2) 
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3. PROTON AND MESON COOLING 

The shock-accelerated protons produce mesons via pp and 
pj interactions. Since the meson multiplicity in pp interac- 
tions is about 1 for pions while 0. 1 for kaons, neutrinos con- 
tributed by pion decay are dominant when the coohng effect 
of pions is not important, which is applicable to the low en- 
ergy pp neutrinos. Therefore we here consider only pion pro- 
duction in pp interactions. The pion production by pp inter- 
action in the sub-photosphere dissipation is efficient since the 
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Fig. 1 . — Inverse of proton cooling time scales in the comoving frame of the 
internal shock as functions of proton energy. The straight solid line, dotted 
line and dashed line are for proton-proton collision, Bethe-Heitler cooUng 
and p7 processes respectively. Also shown is the cooling time scale due 
to adiabatic expansion (which is equal to the dynamic timescale tj^_^^ given 

in the text). The parameters used in the plot are Li^ = lO^^ergs"' , F = 100, 
R = 10"cm and Sj = 300keV. 

coohng time in the shock comoving frame, 

tpp = l/i<Tppn'pCKpp) = 0mSLll,2Ri,rl s, (3) 

can be shorter than the shock dynamic time, t'^^,^ = R/Tc = 
0.03/?iir2'sec, where app = 4 x 10"^^cm^ is the cross 
section for pp interactions, n'p = Lk/{4-KFi^T^mpC^) = 2 x 
IQ^'' Lk.s2R~yi^2^cm~^ is the proton number density, and /Tpp ~ 
0.5 is the inelasticity. Protons also cool through Bethe-Heitler 
interactions (po pe^e') and pj interactions when the tar- 
get photon energy seen by the protons is above the threshold 
energy for each interaction. Denoting by n(e~^)dej the num- 
ber density of photons in the energy range to e^ + dej, the 
coohng time in the shock comoving frame for pj and Bethe- 
Heitler cooling processes are given by 



'{pjJiH} ■ 
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(4) 



where Tp = e'p/mpC^, a and K are respectively the cross sec- 
tion and the inelasticity for p^ (or Bethe-Heitler) process. As 
a rough estimate, the Bethe-Heitler cooling time is 



-7^ =0.5( 



(28/9)hi40-218/27 L^,52^^ e'. 



)(7^)(; 



(28/9)hi2A:-218/27''L^,5i 3keV 



;) (5) 



when k = e'pe'^/inipmec'^) is a large value (a good approx- 
imation when k > 10, Chodorowski et al. 1992), where 
s'^ is the thermal peak energy of photons in the comoving 

frame. So when the proton energy is larger than ^ = 

1500(3^)"'GeV, the Bethe-Heitler cooUng dominates over 
the pp cooling. At even higher energies near the threshold 

for p7 interactions at e^*^' = 6 x 10'*(3|^)"'GeV, p'y coohng 
becomes increasingly dominant. 

We compare these three cooling time scales for protons in 
Fig. 1 for representative parameters, using more accurate cross 
sections for p'y and Bethe-Heitler processes in Eq.(4). For 
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photopion production cross section, we take the Lorentzian 
form for the resonance peak (Miicke et al. 2000) plus a com- 
ponent contributed by multi-pion production at higher ener- 
gies, while for Bethe-Heitler process we use the cross section 
given by Chodorowski et al. (1992). The number density 
of photons used in the calculation has been assumed to have 
a blackbody distribution. The numerical result confirms that 
that pp cooling and p^/ cooling are dominant, respectively, at 
lowest and highest energies, while at the intermediate ener- 
gies, the Bethe-Heitler cooling is the dominant one. 

Defining the total cooling time for protons as t'^ = l/(f'^^ + 
t'gu + t'p^), the total energy loss fraction of protons is rjp = 
Min{f^^,„/f^,, 1}. The fractions of energy loss hy pp and p7 
processes are respectively. 



+ t'~BH + t'~o\) 
-1 



C'7 p-iK^' pp'^^' BH'^^' p-/)- 



(6) 



The cooling of secondary pions may also affect the neu- 
trino production efficiency if they suffer from cooling be- 
fore decaying to secondary products. The pions suf- 
fer from radiative cooling due to both synchrotron emis- 
sion and inverse-Compton emission. The total radia- 
tive cooling time is = 3mlc^ /[4aTm^e'^UBil + fic)] ^ 
0.002(e;/lTeV)-ie^;_iL^i52/;2jr2s, where t/^ is the energy 
density of the magnetic filed in the shock region and fic < 1 
is the correction factor accounting for the inverse-Compton 
loss. The pions also cool due to collisions with protons (Ando 
& Beacom 2005). The cooling for this hadronic process 
is t'^Mad = e'J(c'T.pn;Ae;) = OmeL^'.^Rn^h, where a^p = 
5 X 10~^^cm^ is the cross section for meson-proton collisions 
and Ae^ = 0.5e^ is the energy lost by the meson per collision. 
The suppression of neutrino emission due to cooling of pions 
can be obtained by comparing the cooling time ^.^^ or 

with the lifetime of pions = ^t^t = 1.9 x 10"'*(e^/lTeV)s 
in the shock comoving frame, where '^.^ and r are the pion 
Lorentz factor and proper lifetime. This defines two critical 
energies for pions, above which the effect of radiative cooling 
or hadronic cooling starts to suppress the neutrino flux, i.e. 



<,rad = ^^B^'h%'RnT2T^Y and e'^„„, 
The total cooUng time of pions is f^ ^. = 1 /(C'Lrf + C~L<<^ 
the total suppression factor on the neutrino flux due to pion 
cooling is (Razzaque et al. 2004) 

c. = Min{f;,,/T;,i}. (7) 

4. NEUTRINO FLUX FROM GRBS 

The total energy emitted in neutrinos from pp or p^ pro- 
cesses per GRB is, respectively, 

^.J{pp,n}(^^) - g J^^-^ , (8) 

where Ep is the energy in accelerated protons in one burst 
during the dissipative photosphere phase, e'p ,„„j. and e'p „„„ are 
the maximum and minimum energies of acceleration protons. 
In the absence of pion cooling loss, the neutrinos produced 
by pion decay carry 1/8 of the energy lost by protons to pion 
production, since charged and neutral pions are produced with 
roughly equal probability and muon neutrinos carry roughly 
1 /4 of the pion energy in pion decay ' . The mean pion en- 
ergy is about 20% of the energy of the proton producing the 

' As an approximate estimate, we have neglected tlie effect of multi-pion 
production, the muon contribution decay to the neutrino flux and the neutrino 
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Fig. 2. — Diffuse muon neutrino flux on Earth contributed by pp and p7 
interactions from the dissipative photosphere of GRBs, assuming that energy 
in accelerated protons in one burst is Ep = 1.5 X 10^^^ erg during the dis- 
sipative photosphere phase, the bulk Lorentz factor T = 100, the GRB rate 
i?GRB(0) = IGpc^'yr"' and /; = 3. The solid (red) and dashed (blue) lines 
are for shock dissipation at R = 10"cm and lO'^cm respectively. The dotted 
(green) line is for shock at photosphere radius lO'^cm with a broken power- 
law photon spectrum assumed (see the text for details). The solid lines denote 
the pp neutrino component, while the dashed lines denote the p7 neutrino 
component. Also shown is the Waxman-Bahcall (WB) bound (Waxman & 
Bahcall 1999). 

pion, so the mean energy of neutrinos is ~ Q.Q5ep. As- 
suming that protons are efficiently accelerated in shocks with 
an energy density of Up = lOUi^, the number of TeV neutri- 
nos from one GRB is about Ni, = 0.1($-y/10~'*ergcm"^) for 
77p ~ 1, according to Eq.(8). So only from very strong bursts 
with gamma-ray fluence > 10~^ergcm"^, which are very 
rare events, can neutrinos from single GRB be detected. 

The aggregated muon neutrino flux from all GRBs is ap- 
proximately given by 

'=^l^{pp.pj}('^'^) — (Thk^'^l-^{pp-in}^'^''^^GRBiO)fz 
= 1.5 X 10-%,53(T|^^)(f )?7pC{py.p7}C^GeVcm-2s-isr-', 

(9) 

where fr is the correction factor for the contribution from high 
redshift sources and /?grb(0) is the overall GRB rate at red- 
shift z = 0. Assuming that GRB rate traces the star-formation 
rate in the Universe, the calculation gives f- ~ 3 (Waxman & 
Bahcall 1999). It is not clear how efficiently the protons are 
accelerated in GRB shocks. Assuming an optimistic case that 
protons are efficiently accelerated in shocks and that half of 
the kinetic energy dissipation occurs below the photosphere^, 
we take a mean value Ep = 1.5 x lO^-'ergs for the isotropic 
equivalent energy in accelerated protons in one GRB during 
the dissipative photosphere phase, based on a typically used 
value Lit = lO^^ergs"' for the isotropic kinetic energy luminos- 
ity and a typical long GRB duration of Ar = 30 s. The GRB 
rate^ at redshift z = is taken to /?grb(0) = lGpc"^yr"' (Guetta 

oscillation effect, which may affect the estimate of the pp neutrino flux within 
a factor of 2 (the factor, however, could be larger for p7 neutrino flux). 

^ This is based on the analysis by Ryde (2006) and also in a very recent 
paper by Ryde & Pe'er (2008, arXiv:0811.4135Yl), who find that the thermal 
photons caiTy a fraction of 30% to more than 50% of the prompt emission 
energy. 

' There is lai'ge uncertainty in the estimate of the local GRB rate. Some 
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et al. 2005; Liang et al. 2007). The isotropic luminosity 
is taken to be = 10^'erg The energy-dependent neutrino 
flux contributed by pp and p^ interactions are plotted in Fig. 2 
for a set of representative parameters of the dissipative pho- 
tosphere model and three different dissipation radii. If the ki- 
netic energy is dissipated at radius ofR= 10"cm, the calcula- 
tion (the red solid curves) shows that at energies below tens of 
TeV, the neutrino flux is dominated by a component. Tak- 
ing the detection probability of P^^ = 10"^(e,y/lTeV) for TeV 
neutrinos (Gaisser et al. 1995), the expected flux of upward 
moving muons contributed by this pp component is about 8- 
10 events each year for a km^ neutrino detector, such as Ice- 
cube. We can also estimate the atmospheric neutrino back- 
ground expectation in coincidence with these GRB sources, 
noting that the search for neutrinos accompanying GRBs re- 
quires that the neutrinos are coincident in both direction and 
time with gamma-rays. Taking an average GRB duration of 
~ 30 s, an angular resolution of Icecube of ~ 1 ° and the atmo- 
spheric neutrino background flux of ~ 10~'*GeVcm"^s"'sr"' 
at 1 TeV (Ahrens et al. 2004), the atmospheric neutrino back- 
ground expectation is ~ 5 x 10"^ events from 500 GRBs (in 
one year). Such a low background in coincident with GRBs 
allows the claim of detection of TeV neutrinos from GRB 
sources. At energies from a few TeV to tens of TeV, the 
Bethe-Heilter cooling suppresses the pp cooling, resulting in 
a steepening at several TeV in the neutrino spectrum. At en- 
ergies above the threshold for p^ interactions, the neutrino 
from /77 process is heavily suppressed due to the strong ra- 
diative cooling of secondary pions. For a larger dissipation 
radius at /? = lO'^cm (the blue dashed curves), the neutrino 
emission flux from p7 process is no longer suppressed and 
in this case both pp and p^ neutrino components could be 
detected by km^ detectors. We also calculated the neutrino 
flux, shown by the green dotted lines in Fig. 2, for shock 
at photosphere radius Rph, assuming that the radiation spec- 
trum is a broken power-law (due to Comptonization) peaking 
at = lOOkeV, with lower energy and higher energy photon 
indexes given by -1 and -2 respectively. In this case, the pp 
neutrino flux becomes small (may be marginally detectable), 
while the p^ neutrino flux spectrum is similar to the analytic 
result obtained by Waxman & Bahcall (1997), as expected 
for a broken power-law photon spectrum. Note that in one 
burst the shock dissipation could be continuous from small 
to large radii, as indicated by the larger variability timescales 
seen in GRBs. By comparing the three cases of different dis- 
sipation radii in Fig. 2, one can see that as the shock radius in- 
creases, the neutrino emission from pp component decreases, 
while the p^ component increases until it reaches the satura- 
tion level. The total pp neutrino flux from such continuous 

people suggest a lower GRB rate based on the analysis of Swift bursts with 
Rgrb(O) = 0.05 - 0.27Gpc"Vi"' (e-g- Guetta & Piran 2007; Le & Dermer 
2007) , while others get a higher rate comparable to earlier estimate before 
Swift (e.g. Liang et al. 2007). 

^ Some observations have indicated rather high radiative efficiency, and 



dissipation is thus contributed predominantly by the deepest 
internal shocks below the photosphere. In the whole neutrino 
spectrum, a "valley" is seen between the pp and p^ compo- 
nents of the spectrum, which may be a potential distinguish 
feature of the sub-photosphere dissipation effect. 

5. DISCUSSIONS AND CONCLUSIONS 

Waxman & Bahcall (1997) as well as later works have stud- 
ied the neutrino emission from the photomeson process during 
the prompt internal shocks of GRBs, assuming that the radi- 
ation in the shock region has a broken power-law nonthermal 
spectrum. It was found that the neutrino emission peaks at 
energies above 100 TeV. Towards lower energies, the neutrino 
emission intensity decreases as £^$,y ~ Si, (Waxman & Bah- 
call 1997). However, if internal shocks, especially at the early 
stage of the prompt emission, occur below the photosphere, 
a quasi-thermal spectrum will arise. In this Letter, we have 
discussed the neutrino emission associated with the dissipa- 
tive photosphere that produces such prompt thermal emission. 
We find that pp interaction process becomes important for 
shock-accelerated protons and provides a new neutrino com- 
ponent, which dominates at energies below tens of TeV. The 
neutrino emission from photopion process of protons interact- 
ing with the sub-photosphere radiation could be significantly 
suppressed due to radiative cooling of secondary pions, when 
the dissipation radius is relatively small. Nevertheless, the to- 
tal contribution by photopion process will not be suppressed 
since the shock dissipation could be continuous and occur at 
large radii as well. Although TeV neutrinos may also be pro- 
duced during the early precursor stage of GRB, i.e. before 
the jet breaking out the progenitor star (e.g. Razzaque et al. 
2004; Ando & Beacom 2005), we want to point out that the 
TeV neutrino component discussed here can be distinguished 
from them, because in our case the neutrino emission is asso- 
ciated in time with the prompt emission. 



After this work has been completed and later put onto the 
arXiv website (arXiv:0807.0290), we became aware that K. 
Murase was also working on the sub-photosphere neutrino 
independently (K. Murase, 2008, arXiv:0807.0919). XYW 
would like to thank R Meszaros , S. Razzaque, K. Murase, E. 
Waxman, Z. Li and K. loka for useful comments or discus- 
sions. This work is supported by the National Natural Sci- 
ence Foundation of China under grants 10221001, 10403002 
and 10873009, the National Basic Research Program of 
China (973 program) under grants No. 2007CB8 15404 and 
2009CB824800. 

the importance of pp neutrinos may be reduced if is smaller than lOL-,.. 
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